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Abstract 

We study the renormalizable singlet-doublet fermionic extension of the Standard Model. In this model, the 
new vector-like fermions couple to the gauge bosons and to the Higgs via new Yukawa couplings, that allow 
for nontrivial mixing in the new sector, providing a stable, neutral dark matter candidate. Approximate analytic 
formulae are given for the mass spectrum around the blind spots, where the dark matter candidate coupling to h 
or Z vanishes. We calculate the two particle scattering amplitudes in the model, impose the perturbative unitarity 
constraints and establish bounds on the Yukawa couplings. 
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1 Introduction 


There is a renewed interest in vector-like fermions fT), as these are less constrained than chiral fermions in elec- 
troweak precision tests of the Standard Model (SM) ||2] O and free of anomalies. Vector-like fermions appear 
naturally in various beyond the Standard Model scenarios, such as specific string theories, bulk femions in univer¬ 
sal extra dimensional models 0, little Higgs theories EEl, supersymmetric models as superpartners of standard 
bosons e.g. higgsinos Q , various composite Higgs models Ham and simplified models of dark matter lUoiiniiii. 
Vector-like fermions help gauge coupling unification, but the lower unification scale is in conflict with proton de¬ 
cay constraints Col. A recent survey of phenomenological implications of vector-like extensions HI presents 
constraints from electroweak precision observables, direct collider searches, Higgs production and decays. 

The couplings to the standard light (and heavier) fermions are severely constrained. To avoid flavour problems 
(for flavour issues see m) and to ensure there is no mixing between SM and the new vector-like fermions a Z 2 
matter parity is introduced. The new fermions are odd, while the standard particles are even under this symmetry. 
The lightest new particle is stable and if it is electrically neutral then it provides a dark matter candidate. The 
singlet-doublet dark matter model is an effective theory or simplified version of the neutralino MSSM dark matter 
sector. The mixing in the dark sector allows for a wide range of Higgs and gauge boson couplings that is able 
to avoid direct detection in dark matter experiments and colliders and give the measured relic density, while still 
remaining a relatively simple model. The model contains four new parameters, two dimensionful mass parameters 
and two Yukawa couplings, originally completely unconstrained. A recent analysis in m has found combined 
bounds on the vector masses and the overall value of the Yukawa coupling. 

The spectrum generally is a solution of a third order equation and the analytical treatment is difficult. Another 
problem is that the four-dimensional parameter space can be only visualized if at least two parameters (generally 
the Yukawas) are fixed. In llT4l the phenomenologically allowed regions of the parameter space have been found 
in the neighborhood of the blind spots, where the coupling of the dark matter candidate to h or Z vanishes. In 
these special points one eigenvalue can be found and the third order equation simplifies to a second order one. 
With this observation we have found analytic solution for the masses close to the blind spots. We have further 
shown that the relevant coupling constants are small and can avoid observation by direct detection experiments. 

The parameters of the model can be constrained on theoretical grounds. Perturbative unitarity is a useful tool 
to set limits on effective theories. Starting from weak charged current interaction, it leads to the well-known gauge 
bosons and the Higgs with usual couplings of the SM. It provided an upper bound on its validity without the Higgs 
boson and a theoretical upper bound on the Higgs mass HE). For chiral fermions, there is also an upper bound 
on the scale of fermion mass generation in the few TeV region ifThll . In this paper we calculate the two-particle 
scattering amplitudes involving the new fermions, aiming to constrain the free parameters of the model. As the 
new part of the Lagrangian is renormalizable, the potentially dangerous amplitudes growing with energy cancel 
each other. The elastic scattering of neutral fermions of the doublet and the singlet gives a meaningful amplitude, 
constraining the value of new Yukawa couplings. The new vector-like fermions only partially receive their masses 
from the Higgs, therefore the new bound, <1.23 TeV cannot be directly translated to the mass of the dark 

matter candidate or to the mass parameters of the model. 

In section|2]we review the singlet-doublet vector-like fermion extension of the Standard Model, then comment 
on the current dark matter bounds. In section 0] we give analytic formulae around the experimentally favored 
regions. In section|5]we calculate the new two-particle scattering processes that can contribute to bound the model 
parameters from perturbative unitarity. The paper is closed with conclusion. 


2 The model 


We extend the Standard Model with a pair of SU{2)w doublet Weyl-fermions, tpi = '02 = ( ^0 

that acquire a Dirac mass term together and a singlet, with a Majorana mass term. All three are color-singlet to 
avoid the strong collider bounds, their quantum numbers are listed in Table[T] We also assume a matter parity-like 
Z 2 symmetry that forbids the new fermions to couple directly to the SM ones. The new particles then only couple 
to the Higgs and gauge bosons in pairs, so the Lagrangian can be break up into two parts, C = Csm b- Cds- 
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Table 1: The weak quantum numbers of the new fermions and the Higgs. 


The Lagrangian for the new fermions contains their kinetic and mass terms and Yukawa couplings. 

- ^mdV'iV'2 + ^TOsX°X° + yitpiHx° + 2/2V’2^X° + h.c^ ( 1 ) 

There are two new dimensionful mass parameters and two dimensionless Yukawa couplings, {irid, rUg, yi, 2 / 2 )- 
The three phases of ^ 1^2 and x° can be fixed by setting three parameters (e.g. irid, Ws, yi) to be real and positive. 
We will consider the case where 1/2 is reaQand emphasize that the sign of 2 / 2 / 2/1 is physical. 

A nice and widely used parametrization of the Yukawas is 


yi = ycos9, 7/2 = ysin 6 >. 


( 2 ) 


The decoupled MSSM bino-higgsino system with one light Higgs corresponds to /3 = 0, and y is related to the 
17(1) gauge coupling. 

The charged fermions, and merge into a Dirac fermion, ) with mass nid- Without the 


r2 

Yukawa couplings, the neutral part of the doublets also form a Dirac fermion with mass nid, = 




. As 


the Higgs gets a vacuum expectation value with non-vanishing Yukawa couplings, there is a mixing in the new 
neutral sector. 


with the mass matrix Mn, 


= 


X 




n mAi + h.C; 

(3) 




/ mg 

yiv 

V2 

y2v\ 

V2 \ 


yiv 

V2 

0 

md . 

(4) 

y2V 

\V2 

md 




The corresponding characteristic equation follows. 


{rUg - A) (A^ - ml) + mdyiy2v'^ + A 


(yf + yi) 


= 0. 


(5) 


The cubic equation can be solved analytically using the Cardano formula or numerically. There is a negative 
eigenvalue that can be flipped to be positive, multiplying the corresponding eigenvector by i or equivalently 
performing the Takagi diagonalization on the system and get only the positive masses. Generally the spectrum 
contains three neutral Majorana fermions, the mass eigenstates xi > X 21 X3 ■ The lightest will be denoted by x with 
mass m^, that will be stable due to the Z 2 symmetry. It is an ideal dark matter candidate if lighter than the charged 
fermion < md, made from the following composition of the weak eigenstates, 

X = Uiix° + C/ 12 V'? + \Uii\^ + \Ui2\^ + l«Ti3p = 1. (6) 

Uh characterizes the amount of the singlet in x ■ When md > mg, x is more singlet- or bino-like with >0.5, 
while when md < mg, x is more doublet- or higgsino-like with Uh < 0.5. 

^ It can have a CP-violating phase, which is discussed in Qo). 
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3 Dark matter constraints 


The singlet-doublet model for various values of the parameters can provide a wino- or a higgsino-like dark matter 
candidate. The parameters of the model and the dark matter candidate were recently explored in ifTTl fTSl [191 l20l . 

Several experiments constrain dark matter candidates. Planck observations ED on the relic density abundance 
Q.dmh^ — 0.12. Indirect searches obtained by Fermi-LAT Ii22l and IceCube Ii2^ . and direct dark matter searches 
by PICO EHj lux li25l and XENON 100 In the low mass region, collider bounds are becoming important, 

too, as the charged fermion mass is already excluded for rud < 100 GeV by chargino searches at LEP liZTl and 
the LHC can study the production and annihilation of the dark matter particle li^ or the invisible Higgs and Z 
decays. 

The most stringent bounds are coming from the direct detection experiments in the low mass region. Here 
the relevant couplings of the dark matter are the ones to the Higgs and Z. It is important to note that the Chxx 
coupling is related to the the characteristic equation. After differentiating (|5]l with respect to the Higgs VEV, it 
can be solved for = Chxx- 


^hxx ~ 


{ 2 yiy 2 md + {yf + yl) m^) v 


m‘ 


; + (2/1 + yl) \ 




(7) 


^_ - vl) ~ _ (g) 

2 {m\ - {Ayiy 2 mxmd + {yl + yl) {m\ + m^)) 

where is the eigenvalue, which can be negative as well. 

The spin-independent ’’blind spots” li29l are defined, where Chxx coupling vanishes, similarly the spin-dependent 
blind spots, where czxx coupling vanishes. These points evade the related direct detection bounds. To keep the 
couplings small and avoid overabundance in the relic density, the annihilation cross section needs enhancement. 
This can be achieved, if there is a pole in the s-channel processes, e.g. ttIx ~ ^ or ^ iflTl . A recent work 1(141 
collected all these bounds on the singlet-doublet model, including dark matter searches and colliders. It carries 
out a full numerical study to cover the parameter space, that includes a higher mass region, where vrix > 100 GeV 
and a light mass region, both illustrated with three representative Yukawa couplings y. 

In the large mass region, indirect detection constrains the underabundant region in the parameter space. Eor 
the smaller Yukawas y < O.OI that is the only bound and > 280 GeV with Uii > 0.5 remains available 
for any tan0. Eor the more ’MSSM-like’ y = 0.2, it allows rux ^ 220 GeV with > 0.65. Here the direct 
detection experiments are relevant for thermal x with si 0.12, in the region where iris ~ n^d LUX excludes 
up to uix ^ I TeV mass. In the third region with larger Yukawas y = 1, they obtained uix ^ 275 GeV, unless it 
is purely singlet (Ufi > 0.8). 

In the light mass region x is singlet-like as we mentioned, the allowed regions are around the blind spots with 
nix ~ ^ these are further constrained by the invisible Higgs decay. Parametrically, for the smaller 

Yukawa coupling y = O.I, x is excluded in the 80 GeV < < 220 GeV range with Ui^ < 0.65, and there is no 

further bounds on the y = 1 region. 

The experiments bound the parameter space, but four parameters are challenging to interpret. In the next 
section we explore analytically the regions around the blind spots. 


4 Regions of small couplings 


The bounds from direct detection experiments can be fulfilled with small dark matter couplings to the Higgs 
and Z bosons (|7l[8]l. The LUX li25l and XENONIOO li26ll experiments set bounds on the spin independent (and 
dependent) cross section giving for a dark matter heavier than the nucleon Ii29ll Chxx ^ 0.01 — O.I and similarly 
Czxx — 0 01 — O.I depending on the dark matter mass. In this section we study analytically the region around the 
blind spots, the experimentally favored regions of the parameter space mi. 

The blind spots are the following, 

^Zxx = 0, if \yi \ = \y 2 \ or \nix\ = rud (9) 


and 


Chxx ~ 0, if nix + sm(20)md = 0, 


( 10 ) 
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where sin(20) = ^ ^ . The conditions in (|9l) are not independent. If A = im^ is an eigenvalue of the 

characteristic equation (|5]l it follows that yi — =fj /2 or tan0 = =f 1. Therefore it is enough to study the cases 
in the neighborhood of yi — ±j /2 and (fTOl) . Similarly if Chxx vanishes, it follows from (|5]l that the corresponding 
eigenvalue nix is irrid or uis- In both cases the cubic equation remarkably simplifies by finding one root (m^ or 
nig) and the rest is a quadratic equation. The expansion is performed around these special points. 

Small czxx-> Vi - “ 2 / 2 , tan^ ~ -1 

The mass eigenvalues are 


mi = nid{l — X-), where X- 


jyi +y 2 fv'^ 

4TOd(TOs - md) + ’ 


( 11 ) 


mg-md , 1 /7-^^ 

m 2,3 = -^-± 2 V 


2y2^,2 ^x- I 1 =F 


mg - imd 


xj (Ws+TOd)^ + 2 j/2x,2 


( 12 ) 


where y^ = j/i + j/|. The first two terms in (fT^ are the exact solution of the second order remnant of the 
characteristic equation for nix = md in the czxx — 0 blind spot. The correction is proportional to x-H The 
spectrum in the blind spot j/i = — 2/2 is shown in Fig. [T] 



Figure 1; Mass eigenvalues for yi = — j /2 vs. increasing overall Yukawa terms for nig = 300 GeV < nid = 
500 GeV on the left and nig = 800 GeV > md = 500 GeV on the right. The negative eigenvalue {m 3 ) has been 
flipped to positive value and both are shown. 


For nig < md, the dark matter candidate is the one with mass ni 2 , until y < - \/md{md — mg), then m 2 
becomes larger than mi. Its couplings are 


CtiXX 


Sy'^v 

mg +nid’ 


c^xx 


(2/1 - yDmzv 
2 {ml-ml) ■ 


(13) 


We see that the Z coupling scales according to the blind spot condition (yi + 2 / 2 ), but small Higgs coupling requires 
small overall Yukawa, y and relatively large mg + m^ compared to the Higgs VEV 

For md < mg or y > ^ \/md{md — mg), the dark matter candidate is mostly doublet like and its mass is mi. 
The couplings are 


^ (yi + y2)‘^vmd ^ (yi - yl)mzv 

2md{md - mg) - y^u^ ’ 4md{mg - md) + 2y^v‘^ ' 

Both couplings go to zero with (yi + y 2 )^ and (yi + y 2 ) multiplied with small mass ratios for non-degenerate 
mass parameters \mg — md\ > v. This is the most favored region to avoid direct detection even in the presence of 
sizable Yukawa interactions if they are tuned (yi -f y 2 ~ 0). 

There is a smooth limit for mg = md 


_ iyi+y2)'^md ^ _{yl-yl)mz 

''^XX ~ „,2„. ’ '-■^xx “ 




2 y^v 


( 15 ) 


The Higgs coupling goes to zero with (yi + y 2 )^, while the Z coupling scales with yi + y 2 as before. 


^The real expansion parameter for small Cz^-y in equation is- '^d{vi^V2) " —^^ which has positive definite denominator. 

( 2 m 2 + 6 m 2 + 3 i/ 2 „ 2 )^ 
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Small cz^^, 2/1 ~ 2 / 2 , tan 6 » ~ 1 
The mass eigenvalues are 


mi = —md(l + x+), where x+ = 
rus + md , 1 n -, Wd 


{yi - y 2 fv^ 


4md(ms + rud) + y'^v'^ ’ 

nis + 3md 


m2,3 = -± -\/(m«-md)2 + 2y‘^v'^ + —a;+ 1 T 

II ^ \ ^{nia-rndY + _ 

The spectrum is similar to the first case and it is shown in Fig. |2] 
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Figure 2: Mass eigenvalues for yi = j /2 vs. increasing overall Yukawa terms for nis = 300 GeV < nid = 500 
GeV on the left and ms = 800 GeV > m^ = 500 GeV on the right. The negative eigenvalue (mi) has been 
flipped to positive value and both are shown. 


The lightest fermion with mass m 3 is the dark matter candidate. Its mass starts from ms or m^ and decreases 
as the Yukawa couplings are increased. 

The Z coupling goes to zero with yi — 1/2 in agreement with the blind spot condition, assuming y'^v'^ < 
\m1 — m^l for the sake of simplicity 


and 


~ (yi 2 / 2 ) 


^Zxx ~ ( 2/1 2 / 2 ): 


2 (m 2 - m 2 ) 
mzv 


for ms <md, 


for ms > md- 


2md{ms - md) 

The leading behavior of the Higgs coupling constant for m^ = m 3 is for ms Y md 


(18) 


(19) 


^hxx ~ 


y'^v 


\md -ms\' 


( 20 ) 


The Higgs coupling is small only if the Yukawa mass correction is smaller than the doublet-singlet mass 

splitting. 

For ms = md the couplings are 


, ( 2 /? - 2 /i) mz 

Chxx - 2 / 1 , and ■ 


( 21 ) 


The Z coupling can be small for tuned Yukawas and if the new vectors are heavier than the Z, but the small Higgs 
coupling needs small yi. 


Small Chxxi sin 26* ~ ^ 


As we mentioned the Higgs blind spot is at m^ = imd or m^ = ms- If {m^l — md then | tand| ~ 1 and it is 
discussed in the previous two points. The non-trivial new case is when m^ — ms < md. The mass eigenvalues 
are 


mi = ms(l - z), 


where 


+ y'^v'^ 

2 m^ — 2 ms + 2 /^"*^^ ' 


( 22 ) 
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TO2,3 = ± 


(23) 



The square root is the exact solution of the second order remnant of the characteristic equation in the Chxx = 0 
blind spot for rux = Wg and the correction is proportional to The spectrum is shown in Fig. [3 



Figure 3: Mass eigenvalues for = 0 vs. increasing overall Yukawa terms forms = 300 GeV < rud = 500 
GeV on the left and mg = 800 GeV > trid = 500 GeV on the right. The negative eigenvalue m 3 has been flipped 
to positive value and it coincides with m 2 . 


The coupling constants for the rUx = rui dark matter candidate are 


m 


^hxx ~ ’ ^Zxx — 


{yl - yi)mzv 


Amd {md - rus) + 2y'^i 


(24) 


Small Z boson-dark matter coupling requires additional tuning. The ^ ratio is fixed by therefore for non¬ 
degenerate masses m^ — mg > yv, Czxx scales with In case of rUg = m^ we get back the first 

TTLdy rnl — m^ 

case and both coupling constants (fT 3 T l are small. 

We have seen the exclusions and limits in the parameter space, but there are no bounds on the individual 
parameters. In the following section, we aim to set bounds on the four parameters one by one by exploring the 
consequences of perturbative unitarity. 


5 Perturbative unitarity 

Perturbative unitarity is an essential tool in exploring effective field theories. Amplitudes growing with energy 
indicate the breakdown of the effective theory and the validity range of the model can be estimated or new particles 
and interactions can be added to cancel the terms with bad high energy behavior. In the SM, the complete gauge 
structure of the W, Z and Higgs interactions can be recovered and there are cancellations due to gauge symmetry. 
Unitarity still constrains the parameters of the SM, the Higgs self coupling, which can be translated to the Higgs 
mass 03 . Later, the method was applied to massive chiral fermions without a Higgs boson. It was shown in 
03 that the scattering amplitude of a fermion-antifermion pair to longitudinally polarized gauge bosons must 
be unitarized below 3.5 TeV in the case of the top quark, constraining the scale of fermion mass generation, e.g. 
giving an upper bound on the validity of the model. 

Here, we investigate tree-level elastic two-particle scattering processes. Considering the J = 0 partial-wave 
amplitude and we require perturbative unitarity for a process with scattering amplitude and scattering angle 9. 

1 

ao = -^ J d{cos0)\M\ , iReool < i (25) 

-1 

There are four relevant processes in the model. Charged fermion pair annihilation to W’s, charged and neu¬ 
tral fermion annihilation to W and Higgs, two charged fermions scattering through Z and 7 and finally, neutral 
fermions scattering through Higgs exchange. The amplitudes are included for each helicity channel of the fermions 
for energies much higher than any masses, y/s ^ mw- 


^The expansion parameter for small c^xx s is ^ - '^dVlV2 " ^ ^ 

(2m2+6m2+3;,2«2)5 
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^-(Sl)^+(S2) ^ W-W+ 


The Feynman graphs for the process are shown in Fig. |4] In the t-channel we have three graphs with the three 
neutral mass eigenstates, Xi, 2 , 3 - Since we sum up all possible internal particles, if we do a unitary transformation 
on Xi, 2,3 and sum up all the new states, we get the same result. That means, we can calculate with the electroweak 


eigenstates, ' 1 '° = 




and where couplings are zero. The amplitude is then 




w- 


w+ 


Figure 4: Feynman graphs for T' 4/+ —)> W W~^ scattering. 




o 


V 


^ (^) 

/( 1 -tan" 


In the amplitude matrix, the incoming fermion helicity channels are the following. 


(*i »>)' (+_ +t) 


(26) 


(27) 


As in the case of chiral fermions M, the s-channel amplitudes grow with s, that cancel with the Z and 7 
exchange graphs. The t-channel grows with ^/s, but here it is canceled by the s-channel graphs and there is no 
need for the otherwise absent Higgs exchange. 


^-{s,)Xi{s2)^W-h 

The charged and a neutral fermion can annihilate into a W and Higgs, illustrated in Fig. |5] In the t-channel 
internal lines, we can use again the weak eigenstate Naively using the weak eigenstates as the initial states, 
we find that the W~h scattering process grows with i/i, while there is no problem with the other 

process ^ W~h (the helicities of the fermions are as in (l27l l'). 



Figure 5: Feynman graphs for the ik Xi —kF h scattering. 


Ms,s, ^ W-h) = 


^ + o(^] 


gyiy/s I zr> 

2 mw ^ 


( mw \ 

\V-s ) 


2 mw 


J 


( sin Q 

(^'” 4 '° ^ w-h) = 2/2 


0 


sin 9^ 

2 V 2 


o 


f mw\ 

\V~s) 


(28) 


(29) 
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But taking the mass eigenstates Xi, 2 , 3 , that mixes '1'° and makes the s-channel graph grows with com¬ 
pensating the t-channel. The whole process becomes unitary and leaving the highest order to be proportional to 
5 '^ sin 9u! and a combination of the mixing matrix elements. 

^-(Si)^+(S2) ^ ^-(S3)^+(S4) 

The charged fermions can scatter through s- and t-channel Z and 7 exchange illustrated with Feynman graphs in 
Fig.| 6 ] 





/ 9 ^ sin^ I 

2 cos^ 

0 

0 

cos^ I 
\ 2 cos^ 9^ 


0 

0 

-g^ 

2 cos^ 

0 


0 

2 cos^ Su, 
0 

0 


g^ cos^ I \ 
2 cos^ 9u 
0 


2 cos^ 



(30) 


For the two-to-two fermion scattering amplitudes the sixteen fermion helicity channels are in the following order. 


/ - 


(si 


S 2 S3 



- + -+\ 



(31) 


In the first two processes, the growing terms in the amplitude cancel each other because of the gauge symmetry. 
The remaining term, just as in the previous two processes, is proportional to the gauge coupling and satisfy the 
unitarity constraint (l25T l. 


^°(S1)X°(S2) ^ ^°(S3)X°(S4) 

The neutral fermion scattering can be still interesting, since this process includes Higgs boson exchange and the 
corresponding unconstrained Yukawa couplings. The related Feynman graph is shown in Fig. |7] This process 
is basically the XiXj XkXi mass eigenstate scattering. But since there is only one type of graph contributing, 
which is proportional to some combination of the Yukawa couplings, 2 , we can mix these processes to get the 
vp^x^ weak eigenstate scattering with the following amplitude. 



Figure 7: Feynman graphs for ^P^x^ —>■ T'°x° scattering. 
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iM («'°,bi)x°2(P2) ^ h{q,) ^ 'i'°s{P3)x°AP^)) = 

■ \2 

1 \ 1 

(y+ + y-75)Msibi))-2 (“safe) iy+ + y-X5)vsAP4)) (32) 

2V2/ s-rrif^ 

Where the shorthand notations of the couplings are y± = j/i ± ?/2 ■ 

0 0 0 \ 

0 0 Oy 

The fermion helicity channels are in the same order as in the charged fermion scattering in (ISTT i. 

Now using oo partial-wave unitarity (l25T l for the non-zero elements in the amplitude matrix, we get a bound 
on the Yukawa couplings. 

|yi,2| < 4v^ Ri 7.1 |yit/2 | < IStt ft! 50.3 (34) 

The Yukawa contributions to the fermion masses are bounded by 

< 2V^v ft 1.23 TeV, (35) 

V2 

where v = 246 TeV is the Higgs VEV. 

We see that perturbative unitarity sets meaningful bounds on the so far unconstrained new Yukawa couplings. 
The other processes are proportional to the gauge coupling and the bounds are automatically satisfied. 

6 Conclusion 

We have studied the vector-like fermionic singlet-doublet extension of the Standard Model. This is a minimal, 
gauge invariant and renormalizable model, motivated by dark matter. It has four free parameters, two masses and 
two Yukawa couplings. A Z 2 matter parity is introduced, the lightest neutral mass eigenstate can be the dark 
matter candidate. As a result of the mixing between the doublet and the singlet there are tree-level couplings of 
the dark matter to the Z and Higgs bosons and the couplings can vary in a wide range. This model is a consistent 
simplified version of UV complete supersymmetric theories, such as the bino-higgsino sector of the MSSM or the 
singlino-higgsino sector of the naturalness motivated NMSSM-like models. The main difference from the MSSM 
bino-higgsino system is that the Yukawa couplings are free parameters here and not related to the the hypercharge 
coupling and tan /3. The effective or simplified model is ideal to test the parameters with experiments. 

To have a better view and analytic control of the parameters, we have calculated the mass eigenvalues and the 
relevant Z, h couplings in the neighborhood of the blind spots. We identified the best region allowed by direct 
detection experiments, where rud < rUg and the not necessary small Yukawas are tuned to nearly cancel each 
other 2/1 -f 2/2 — 0. There are other blind spot regions, where the Yukawas {y) must be small. Direct and indirect 
dark matter searches constrain the combined parameter space, but not the individual parameters of the model, 
still leaving room for a ’’WIMP miracle”. To set bounds on separate parameters, we have calculated two particle 
scattering amplitudes in the model. Applied the bounds of perturbative unitarity and found that no amplitude has 
bad high energy behavior as expected from renormalizability. There are no constraints on the Dirac- and Majorana- 
mass. The new Yukawa couplings appear in s-channel Higgs exchange graphs and are bounded \yi ^2 \ <7.1. The 
LHC phenomenology of the model is studied e.g. in lllllll. The medium and the low mass region can be tested 
at the i/s = 13 TeV LHC in the next few years llT4l . 


Mg.g^ggg, (4/V ^ 'I/V) = 


/O 

0 

0 

Vo 
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